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We show that a topological phase supporting Majorana fermions can form in a 2D electron gas 
(2DEG) adjacent to an interdigitated superconductor-ferromagnet structure. An advantage of this 
setup is that the 2DEG can induce the required Zeeman splitting and superconductivity from a single 
interface, allowing one to utilize a wide class of 2DEGs including the surface states of bulk InAs. 
We demonstrate that the interdigitated device supports a robust topological phase when the finger 
spacing A is smaller than half of the Fermi wavelength Af- In this regime, the electrons effectively 
see a "smeared" Zeeman splitting and pairing field despite the interdigitation. The topological phase 
survives even in the opposite limit A > Af/2, though with a reduced bulk gap. We describe how 
to electrically generate a vortex in this setup to trap a Majorana mode, and predict an anomalous 
Fraunhofer pattern that provides a sharp signature of chiral Majorana edge states. 



Topological superconductors have attracted consider- 
able recent interest because they may provide the first 
unambiguous realization of Majorana fermions in any 
physical setting. The pursuit of these elusive objects 
in condensed matter [1, 2] is motivated largely by the 
non-Abelian statistics [3-5] that they underpin, which 
is widely sought for quantum computation [6] . Although 
much attention recently has focused on finding Majorana 
fermions in ID systems [7-12], 2D platforms [3, 13-18] 
offer some unique virtues such as the ability to perform 
interferometry [19-25] to probe non-Abelian statistics. 
One promising 2D scheme involves a quantum well sand- 
wiched between an s-wave superconductor and a mag- 
netic insulator [15]. Fabricating this device is, however, 
rather nontrivial as one must synthesize high quality in- 
terfaces on both sides of the quantum well — which is typ- 
ically buried in a heterostructure. One can avoid a mul- 
tilayered architecture by invoking a specific type of 2D 
electron gas (2DEG) with appreciable Rashba and Dres- 
selhaus coupling [17], but the candidate materials for this 
proposal are limited. 

In this manuscript we introduce a new 2D Majo- 
rana platform (Fig. 1) consisting of interdigitated su- 
perconductor /ferromagnet insulator strips deposited on 
a 2DEG (periodically modulated ID topological super- 
conductors were considered in Ref. [26]). The pro- 
posed setup exhibits several virtues. For one, our de- 
vice requires interface engineering on only one side of 
the 2DEG — alleviating one experimental challenge with 
previous semiconductor-based proposals. Because of this 
feature one can also employ a wider variety of 2DEGs, 
including surface states of bulk semiconductors such as 
InAs [27-29]. Meanwhile, this structure naturally al- 
lows one to electrically generate vortices to trap Majo- 
rana zero- modes, potentially allowing Majorana fermions 
to be braided using currents similar to the proposal of 




FIG. 1. Schematic of the proposed interdigitated 
superconductor-ferromagnet 2DEG architecture. The device 
supports a robust topological phase when the finger spacing A 
is smaller than half of the Fermi wavelength (i.e. X < Af/2). 
A topological phase can also appear in the regime A > Af/2, 
though generally with a suppressed gap. Table I provides val- 
ues of Af for select 2DEGs when the chemical potential is set 
to fi = 0. 



2DEG 


a [eVA] 


m/rrie 


Af [M m ] 


InGaAs/InAlAs [31] 


0.05 


0.04 


1.19 


InSb/InAlSb [32] 


0.14 


0.0139 


1.22 


Bulk InAs surface [27-29] 


0.11 


0.03 


0.72 



TABLE I. Effective mass m in units of the electron mass ra e , 
Rashba coupling strength a, and Fermi wavelength Af eval- 
uated at \i — for the 2DEG's listed in the left column. 



Ref. [30]. We further show that in our device (as well 
as any 2D topological superconductor) Majorana edge 
states can be detected by observing an anomalous shift 
of the zeros in the Fraunhofer pattern measured in a long 
Josephson junction. 

We model the semiconductor in this device with the 
following Hamiltonian, 
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where creates an electron with spin a and effective 
mass m, /i is the chemical potential, a is the Rashba 
coupling strength, and a a are Pauli matrices that con- 
tract with the spin indices. The spatially varying Zee- 
man and pairing fields induced by the alternating fer- 
romagnetic and superconducting strips are respectively 
denoted by V z (r) and A(r). For simplicity we will re- 
tain only their maximal Fourier components and take 
V z (r) = 2V z s'm 2 (lQx) and A(r) = 2Acos 2 (^Qx). Here 
V z and A are the spatial average of these quantities, 
which modulate at wavevector Q — 2ir/X with A the fin- 
ger spacing shown in Fig. 1. This choice is expected to 
not only quantitatively capture the effects of interdigita- 
tion, but as we will see also leads to an intuitive physical 
picture for the device's behavior. 

As a primer it is worth recalling the physics of the 
sandwich structure originally proposed by Sau et al [15], 
where a uniform Zeeman field V^ umf opens a chemical po- 
tential window in which only one Fermi surface is present. 
Incorporating s-wave pairing with strength A unif in this 
regime effectively drives the 2DEG into a topological 
p + ip superconductor due to the interplay with spin- 
orbit coupling. [2, 15, 17] Quantitatively, the topological 
phase appears provided (V^ unif ) 2 > (A unif ) 2 + \i 2 . In 
our interdigitated setup it is natural to expect that when 
the Fermi wavelength Xp for the semiconductor greatly 
exceeds the finger spacing A, electrons in the 2DEG ef- 
fectively experience 'smeared' Zeeman and pairing fields 
with strength V z and A. Similar physics to the uniform 
case ought to then emerge — in particular, a topological 
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phase when V z > A + fi 2 . 

To confirm this intuition and extract the phase di- 
agram for arbitrary Xp/X, we study the quasiparticle 
spectrum for Eq. (1). Defining a Nambu spinor \I/k = 
[^(k), ^(k), ^|(— k), — k)] T , the Hamiltonian can 
be written in momentum space as 
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with r a Pauli matrices that act in particle-hole space. 
The Hamiltonian describes a semiconductor proxi- 
mate to a uniform superconductor and ferromagnet and 
is precisely the model studied in Ref. [15]. The bulk ex- 
citation spectrum obtained from in the topological 
phase with /i = 0, V z = 1.5 A, ma 2 = 3 A and k y = 
appears in the red dashed lines of Fig. 2; roughly, the 
gap at k x = is set by V z while A determines the gap at 
the Fermi wavevector kp = 2i\jXp. Our interdigitated 




FIG. 2. Bulk quasiparticle spectrum versus k x in a uni- 
form structure (dashed curve) and interdigitated device with 
Af/A = 2.6 (solid curve). In both cases we use parame- 
ters k y = 0, n = 0, Vz = 1.5A, and ma 2 3A. The 
red arrow indicates the pairing gap at the Fermi momentum 
kF = ~ 2ma/h 2 , while the blue arrow denotes the degen- 
eracy gap opened at k x = j due to the interdigitation. Note 
that the excitation spectra for the uniform and interdigitated 
systems differ appreciably only at rather higher energies here. 



structure produces a new term SH^ that couples spinors 
with wavevectors k and k ± Qx. As we 'turn on' these 
couplings the spectrum of evolves very similarly to 
the band structure of free electrons in a weak periodic 
potential [33]. In particular, the dominant effect of 
is to open a gap in the excitation spectrum whenever 
the energies cross Bragg planes at k x = ±Q/2 = ±7r/A 
(modulo reciprocal lattice vectors). For momenta away 
from these values SH^ couples states that are far from 
resonant and hence perturbs these only weakly. 

It follows that for Xp/X ^> 1 the periodic modulation 
modifies the quasiparticle spectrum appreciably only at 
very high energies. This point is illustrated by the solid 
curves in Fig. 2, which display the numerically obtained 
spectrum for the full Hamiltonian in Eq. (2) in a repeated 
zone scheme, using the same parameters as above but 
now with Xf/X = 2.6. Even for this ratio of Xp/X, the 
spectrum is nearly identical to that of the uniform case 
away from k x = ±7r/A. When Xp/X ^> 1 one can clearly 
incorporate 5Hk while essentially leaving the bulk excita- 
tion gap exhibited by the uniform system intact. Thus by 
adiabatic continuity our interdigitated device supports a 
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topological phase in this limit provided V z > A + /i 2 , 
consistent with the intuition provided earlier. As fur- 
ther evidence, Figs. 3(a) and (b) display the quasiparti- 
cle spectrum as a function of k y in a system with open 
boundary conditions along the x direction. The data 
correspond to \i = 0, ma 2 = 3A and Xp/X = 2.5, while 
the Zeeman energy changes from V z = 0.5 A in (a) to 
V z = 2 A in (b). In (a) a trivial gapped state clearly 
emerges due to the weak Zeeman energy. The larger V z 
value in (b), however, satisfies our topological criterion, 
and one indeed sees the signature gapless chiral Majorana 
edge states inside of the bulk gap. 

As one reduces the ratio Xp/X to a value of order one 
or smaller, the physics becomes considerably more subtle. 
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FIG. 3. Quasiparticle spectrum in various regimes for an in- 
terdigitated device with periodic boundary conditions along 
y but open boundary conditions along x. In all parts we take 
\i — and ma 2 = 3 A, while the finger spacing and Zeeman 
energy vary as (a) Af/A = 2.5, V z = 0.5A, (b) Af/A = 2.5, 
T4 2 A, and (c) Af/A = 1.5, T4 = 2 A. A trivial state 
appears in (a) while the larger Zeeman field in (b) drives 
a topological phase supporting chiral Majorana edge states 
within the bulk gap. Interestingly, the topological phase and 
associated Majorana edge states survive even in (c) despite 
the relatively small ratio of Af/A. 



Indeed, once Af/A « 2 the Bragg plane at k x =Q/2 ap- 
proaches the Fermi wave vector, and the pairing gap can 
then be dramatically altered by the interdigitation. We 
ascertain the global phase diagram of our device by nu- 
merically computing the minimum excitation gap S for a 
system on a torus as a function of Af/A and V z /A. Fig- 
ure 4 shows the results for /i = and spin-orbit energies 
of ma 2 = 3. 2 A in (a) and ma 2 = 1.3 A in (b). The follow- 
ing points are noteworthy here: 1) At 'large' Af/A topo- 
logical superconductivity appears when V z > A, in line 
with our results above. 2) The topological phase survives 
over a range of parameters even for rather small values of 
Af/A, though the gap is generally reduced compared to 
the large Af/A limit. Figure 3(c) illustrates the spectrum 
in the Af/A < 2 regime for a system with open boundary 
conditions along x; just as in Fig. 3(b) the characteristic 
chiral edge states again appear here. 3) Interestingly, for 
Af/A ~ 1 the critical value of V z required to generate 
the topological phase decreases compared to the uniform 
case. 

Having numerically demonstrated that our device ex- 
hibits a topological phase with an edge state, we now 
describe how the interdigitated structure naturally al- 
lows us to electrically generate vortices to trap Majorana 
zero- modes. Consider the setup of Fig. 5. Supercurrent 
flowing from contact 1 to contact 2 produces a winding 
in the superconducting phase 0(r) across the fingers in 
the device. When the phase difference between the con- 
tacts approaches 27r, a vortex forms near the center of 
the system to minimize the energy E oc J d 2 r(\70) 2 . The 
spatial profile of the phase follows from the supercur- 
rent j(r) oc A*(r)VA(r) - A(r)VA*(r), with A(r) = 
Asc{ Y )^ ie<KV \ where the pairing potential's magnitude 



( a ) ma 2 3.2 A s/A ^ rna 2 = 1.3 A 8/ A 




FIG. 4. Phase diagrams for \i — and spin-orbit energies (a) 
ma 2 = 3. 2 A and (b) ma 2 = 1.3 A. The horizontal axis repre- 
sents ratio of the Fermi wavelength Af to the finger spacing 
A, while the vertical axis is the Zeeman energy normalized by 
pairing strength. The shading indicates the bulk gap S nor- 
malized by A. Red dashed lines denote the boundary between 
topological phase and trivial phases. 
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FIG. 5. (a) Scheme to electrically stabilize a vortex binding a 
Majorana zero-mode. Here the singular phase winding is in- 
duced by current flowing from contact 1 to contact 2, rather 
than from a magnetic field. In (b) we illustrate the probability 
density extracted from the near-zero energy mode generated 
by a current-induced vortex at the center of the device (pa- 
rameters are /x = 0, ma 2 — 1.3 A, V z — 2 A, and A — Af/4). 
The large central peak corresponds to the Majorana bound to 
the vortex, which hybridizes weakly with the outer Majorana 
running along the perimeter. 

satisfies Asc( T ) — 2 A beneath the superconductors [blue 
regions in Fig. 5(a)] and goes to zero under the ferromag- 
nets [pink regions in Fig. 5(a)]. In particular, one can ex- 
tract 0(r) by iterating the current conservation equation 
V • j(r) = subject to boundary conditions along the 
system's perimeter. With this phase in hand, one can 
diagonalize the Hamiltonian in the presence of a current- 
induced vortex and extract the wavefunctions for each 
quasiparticle state. Figure 5(b) illustrates the resulting 
probability distribution for the near-zero-energy state in 
the spectrum; the large central peak corresponds to a lo- 
calized Majorana mode bound to the vortex core, while 
the outer peak represents a second Majorana mode run- 
ning along the edge. 

Finally, we discuss the detection of Majorana edge 
states in the topological phase exhibited by our device 
(or, equivalently, any other realization) via an unconven- 
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tional Fraunhofer pattern. [34] Consider a pair of topo- 
logical superconductors forming a long Josephson junc- 
tion of width w pierced by a magnetic field [see Fig. 6(a)]. 
At low energies it suffices to focus only on the chiral edge 
states, which can be modeled by an effective Hamiltonian 
H = H t + Hi) + i^tunneiing- [3, 25, 35] The first two terms 
H t / b = ±ivh j dxj t /bd x j t /b describe the kinetic energy 
for the top/bottom edge states, with j t / b Majorana op- 
erators and v the edge velocity. The last term incorpo- 
rates inter-edge tunneling with strength t at the interface 
and reads ^tunneling = it dxjtjb cos[0(a:)/2], where 

6{x) is the local superconducting phase difference across 
the junction induced by the applied field. Neglecting the 
magnetic field that is produced from the tunneling cur- 
rent, Q(x) is determined by the external magnetic flux $ 
according to 6(x) = Oq + 2-7T J^- ^ ($o is the flux quantum 
and #o is the phase difference at the junction's center). 

The Majorana- mediated contribution to the local cur- 
rent density flowing across the junction follows from 
j(x) = jr- sm[0(x)/2]i / y t j b . We calculate the current 
perturbatively in t assuming the weak-tunneling limit 
< 1 where the hybridization energy is smaller than 
the level spacing. In this case the physics depends sharply 
on whether, at t = 0, Majorana zero- modes exist at each 
edge. If neither edge supports a zero-mode, then the cur- 
rent vanishes to first order in t. However, if zero-modes 
exist at both edges (due to an odd number of vortices in 
their bulk) then a finite current (j(x)) = j^sm[0(x)/2] 
emerges, where L is the superconductors' perimeter. In- 
tegrating over the junction width yields a total Major ana- 
mediated current 
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The solid black curve in Fig. 6(b) illustrates \Im\ as a 
function of remarkably, the zeros occur at even multi- 
ples of ^/^o in contrast to the conventional Fraunhofer 
pattern shown for comparison in the red dashed curve. 
For a sample of size 5 jam x 5 jam with the coupling en- 
ergy t = 0.025meV [25], we estimate that the typical 
magnitude of Im is ~ 1.5nA. This result is valid when 
the edge velocity obeys v > 3 x 10 4 m/s so that the weak 
tunneling limit is satisfied. It is important to keep in 
mind, however, that the experimentally observed current 
will not be given by Im alone — a potentially much larger 
conventional current I s flows in parallel. The magnitude 
of the total current I t ot = Is J r Im is sketched by the blue 
curve in Fig. 6(b). One can infer the existence of Im 
by the unconventional Fraunhofer pattern that exhibits 
shifted zeros as shown in the figure. We note that very 
recently an experiment of this type has been performed 
in a long Josephson junction formed at the surface of 
a 3D topological insulator [36], though the findings are 
rather different from what we predict here. 

In conclusion, we have shown that our interdigitated 



Topological 
superconductor 
— — ■ > ■ ■ 



(b) 



»7t 
B 



l< I I L 



Topological 
superconductor 
> 



IP 

If 

HZ. 



lb 



M 

ll 



-10 1 
*/*0 



FIG. 6. (a) Long Josephson junction formed by adjacent topo- 
logical superconductors. A magnetic field B orients perpen- 
dicular to the plane and uniformly penetrates through the 
junction, (b) The solid black curve represents the magnitude 
of the tunneling current arising from coupled Majorana zero- 
modes at the edge. This contribution exhibits zeros at even 
multiples of the flux quantum in sharp contrast to the Fraun- 
hofer pattern exhibited by ordinary s-wave superconductor 
junctions (red dashed curve). The blue curve represents the 
anomalous Fraunhofer pattern that would arise in an experi- 
ment due to the Majorana-mediated component and a parallel 
conventional current contribution. 



structure exhibits a topological phase that is particularly 
robust when the finger spacing is smaller than half of the 
Fermi wavelength. There the bulk gap can be comparable 
to that in a uniform system; furthermore, additional per- 
turbations induced by the interdigitation (such as varia- 
tions in chemical potential and Rashba strength) should 
play a minor role. To access this regime the finger spacing 
should be < 600nm for the quantum wells listed in Table 
I and < 400nm for the surface state of bulk In As. We also 
note that since electrons effectively see 'smeared' fields in 
this limit, the specific interdigitated pattern studied here 
is by no means required — similar physics should arise, 
e.g., in checkerboard arrangements. An interesting fea- 
ture of our setup is that vortices can be generated by ap- 
plying currents. This mechanism may eventually provide 
a practical means of manipulating and braiding vortices 
for quantum computation. We also pointed out that chi- 
ral Majorana edge states produce an anomalous Fraun- 
hofer pattern that can be observed in any realization of 
topological p + ip superconductivity. 
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